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ABSTRACT 

Eclipses in the 14.4 day period double-lined binary SS Lac were observed 
photographically and visually early in the 20 th century, but stopped some 50 or 
60 years ago. This has been explained by the presence of a distant third star 
in the system, which has now been detected spectroscopically with a period of 
679 days. The plane of the orbit of the binary is changing relative to the line 
of sight in response to perturbations from this third object. A recent analysis 
by Milone et al. of all photometric material available for the system, including 
a re-measurement of original Harvard plates, has confirmed earlier reports of 
changes in the depth of the eclipses as a function of time, which are due to 
the third star. In this paper we discuss our detailed analysis of the eclipse 
amplitude measurements, and extract from them information on the change in 
the inclination angle of the binary over the last century. Our use of a much 
improved ephemeris for the system by Torres & Stefanik was found to be crucial, 
and prompted us to re-determine all the amplitudes from the historical data 
at our disposal, including the Harvard material used by MOO. Systematically 
lower measurements on the branches of the minima were properly accounted 
for, and we made use of both a linear approximation to the time variation 
of the inclination angle and a more realistic model based on the theory of 
three-body interactions ("regression of the nodes" effect). The nodal cycle is 
found to be ~600 yr, within which two eclipse "seasons" occur, each lasting 
about 100 yr. The non-eclipsing status of the system is expected to continue 
until the beginning of the 23 rd century. 

Subject headings: binaries: eclipsing - - binaries: spectroscopic — stars: 
fundamental parameters — stars: individual (SS Lacertae) 
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1. Introduction 

The cessation of eclipses in a binary star is a rare phenomenon that can most often 
be explained by the presence of a third object in the system inducing perturbations in the 
orbital elements of the inner pair. The modulation of the inclination angle to the line of 
sight, in particular, is the direct cause of a change in the orientation of the plane of the 
binary such that eventually the stars no longer block each other's light at conjunction. Few 
examples of this phenomenon are known, among them AY Mus (Soderhjelm 1974) and 
V907 Sco (Lacy, Helt & Vaz 1999). 

An especially interesting case is SS Lacertae (HIP 108981, BD+45°3782, P = UA days, 
A3 V, a = 22 /l 04 m 4r!6, 5 = +46°25'38", epoch and equinox J2000), which stopped eclipsing 
around the middle of the 20 th century. Our knowledge about light variations in this binary 
when it was eclipsing comes entirely from visual and photographic measurements that go 
back more than a century, but that are unfortunately of rather poor quality. However, 
the object has received considerable attention recently as a result of a re-measurement of 
archival plate material and new spectroscopic observations. Lehmann (1991) re-measured 
original plate material from the Sonneberg Observatory (1890-1989), presenting the first 
evidence that the depth of the eclipses had changed over the years. He correctly interpreted 
this as due to the presence of an unseen third star in the system, which is gradually changing 
the inclination angle of the inner pair. Tomasella & Munari (1998; hereafter TM98) derived 
a double-lined spectroscopic orbit for the binary, confirming earlier reports that the system 
had not been disrupted by a chance encounter with another star and that the period had 
not changed, as had been proposed initially to explain the disappearance of eclipses. Milone 
et al. (2000; hereafter MOO) carried out an exhaustive analysis of all available light curves 
of SS Lac, and presented clear confirmatory evidence of changes in the depth of the eclipses 
using their new measurements of the original Harvard Observatory patrol camera plates (see 
Dugan & Wright 1935; hereafter DW35) as well as other published measurements. Further 
extensive spectroscopic observations were presented by Torres & Stefanik (2000; hereafter 
TS00), who showed that the system is indeed triple. They determined the orbital elements 
of the distant third object, with a period of about 679 days and a slightly eccentric orbit. 
They also detected apsidal motion for the first time, which is another manifestation of 
perturbations caused by the third star. Further properties and the history of observations 
of the system are described in detail by MOO and TS00. 

Adding to its interest, SS Lac is a member of the open cluster NGC 7209, and therefore 
other information such as estimates of the age, distance, and metal abundance is available. 
Conversely, accurate dimensions for the binary components, if they can be determined, offer 
a chance to perform critical tests of stellar evolution models as well as valuable clues on the 



-4- 



cluster itself, since SS Lac is located very close to the turnoff point in the H-R diagram. As 
pointed out by TSOO, their high-quality spectroscopic orbital solution for SS Lac has the 
potential of giving very precise absolute mass determinations if the inclination angle at the 
current epoch can be established accurately enough. Spectroscopic observations only allow 
the minimum mass of each star to be determined, in the form of the quantity Msin 3 i. The 
value of the inclination at the present time can only be inferred by extrapolation, and this 
requires in particular that its rate of change be known very accurately. Estimates have 
been made in a variety of ways by a number of authors, as described later, but are typically 
based on very limited data. 

The motivation for this study is precisely to use the new information that has now 
become available in the form of light curves analyses and measurements of the amplitude 
of the eclipses, to improve our knowledge of the time dependence of the inclination angle 
and thus to allow the absolute masses and also the radii of the component stars to be 
established. We explore ways of doing this and apply a physically realistic model based on 
the known properties of the binary and the effects of the third star expected from theory. 
Although as we describe later this potential can still not be realized because of limitations 
in the observations, new insights can be gained on the timescale and the pattern of the 
appearance and disappearance of eclipses. 

2. Modeling the change in the amplitude of the eclipses 

MOO presented measurements of the depth, D, of the eclipses in SS Lac that show 
convincingly that the amplitude of the minima decreased gradually from the 1890 's and 
early 1900's to shortly before the middle of the 20 th century, when eclipses were last 
recorded. The peak amplitude is believed to have occurred around 1910 or so. We 
will re-examine these data in detail in §3, but we focus first on a procedure to extract 
information on the rate of change of the inclination angle from these measurements. 

As a means of estimating the approximate time of onset and cessation of eclipses, or 
the beginning and end of the eclipse season, as they termed it, MOO presented a quadratic 
fit to the measurements around the primary eclipse as a function of time (their Fig. 7) that 
provides a reasonable representation of the data. They correctly pointed out, however, that 
a quadratic fit is not necessarily a realistic model for the true variation, which is more likely 
to taper off at either side of the maximum rather than end abruptly. As an alternative, 
MOO discussed a fit to a model involving the complementary error functionF], which has the 



Due to a misprint the equation appearing in MOO is not the actual model they used, which should read 
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desired property of decreasing slowly at both ends although it still lacks any physical basis 
and is merely a mathematical tool. 

The true variation of the depth of the eclipses may be separated into two parts: a 
function we call D(i) (one for each minimum), which describes the change in the amplitude 
as a function of the inclination angle i, and a second function i(t) that describes the time 
variation of %. This latter function is what we seek to determine from the observations. 
We will assume here that the inclination angle is the only orbital element that changes 
significantly with time due to the third star, although strictly speaking this may not be 
true for SS Lac. On the one hand apsidal motion has been detected (TSOO), and the secular 
variation in u (the longitude of periastron) will alter the shape of the light curve slightly 
over long timescales (~1000 yr; see TSOO). On the other hand the perturbations from the 
third star can also cause a modulation of the eccentricity of the binary (see, e.g., Mazeh 
& Shaham 1979), although this effect has not yet been detected. Neither of these effects 
will change the amplitude of the eclipses appreciably in this case, so for our purposes all 
elements except i may be assumed to be constant. 

Similar arguments to describe the time variation of D appear to have been used by 
Lehmann (1991) in the first report presenting evidence of morphological changes in the 
light curve of SS Lac, although few details are given. He correctly attributed the changes in 
the inclination angle to perturbations from a third object in the system, and estimated the 
timescale for the variation. 

2.1. Determining D(i) 

The change in the depth of the eclipses as a function of the inclination angle can be 
derived in a fairly straightforward way once a fit to a representative set of observations has 
been made, using any of the current computer codes to analyze a binary light curve. One 
simply holds all elements fixed and changes % to produce a family of synthetic light curves, 
and then measures the depth of the eclipses (separately at both minima) on each of these 
as a function of i. A table can then be produced that represents D(i), and this may be 
interpolated for ease of use. The critical step here is arriving at a light curve solution from 
the real observations that accurately represents the physical properties of the components. 



y = a ■ crfc { \{x — 6)/c] 2 }. The quadratic argument makes the function symmetrical around the maximum. 
For the benefit of other readers we note also that the phases listed in their Table 10 and Table 11, while 
computed from the TM98 ephemcris, have an integer number of cycles added (2713) to render them positive. 
Also, the cycle number for their primary minimum =^40 should read 2421.988. 
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The very fact that the light curve is changing with time and the relatively poor quality 
of the historical brightness measurements of SS Lac make this step more difficult than it 
would seem. 

Light curve analyses of the photometric observations for SS Lac have been described in 
detail by MOO. They used the University of Calgary version of the Wilson-Devinney code 
(Wilson 1992; Milone, Stagg & Kurucz 1992; Stagg & Milone 1993; Kallrath et al. 1998), 
and examined the three most complete data sets available: 

1. Photographic measurements by DW35 based on Harvard College Observatory patrol 
camera plates, covering the interval 1890-1934. These were published only in graphical 
form as an average light curve, already folded using the ephemeris by the authors. A 
sizeable fraction of these plates and subsequent plates taken at Harvard were re-measured 
by Schiller, Bridges, & Clifton (1996), although these new measures were not used by MOO 
for their light curve analyses. 

2. Photographic measurements by Wachmann (1936) based on plates from Bergedorf 
spanning the interval 1924-1932. Low light levels occur only during the last two years of 
this interval. 

3. Visual estimates by Kordylewski, Pagaczewski, & Szafraniec (1961) between 1927 
and 1948. 

In all cases the solutions were iterated with fits to the radial velocity measurements 
of TM98, under various assumptions on the relative temperatures, luminosities, and mass 
ratio of the components. In their attempt to avoid smearing of the light curve due to real 
changes in shape as a function of time, MOO analyzed these three data sets separately. As a 
result some of them are rather sparse, which adds to the difficulties imposed by the poor 
quality of the observations. In addition, there are some inconsistencies between separate 
solutions, and in the end MOO adopted their results for the DW35 data set as being more 
reliable. 

A light curve solution for SS Lac was also presented by TS00, who chose to merge 
the datasets by DW35 and Wachmann (1936) since they overlap in time. They used the 
computer program EBOP (Etzel 1981; Popper & Etzel 1981), which, although simpler than 
the Wilson-Devinney code, is quite adequate for this well-detached system. 

The MOO and TS00 solutions show some significant differences, no doubt a reflection 
of the difficulties mentioned above, but also because the MOO analyses were done without 
the benefit of the new spectroscopic information obtained by TS00, which only appeared 
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later. The light ratio (secondary /primary)^ derived by MOO is approximately 2.2, while the 
spectroscopic measurements by TSOO indicate a value quite close to unity (1.06 ± 0.03), 
with much larger ratios being clearly ruled out by those observations. TM98 seem to concur 
with TSOO, at least qualitatively, reporting that the secondary is "slightly" more luminous 
than the primary from the relative depths of the Doppler cores and Ha line wings in their 
spectra. Because of the inherent difficulty of determining the ratio of the radii, k, from 
light curves such as that of SS Lac in which the stars are nearly equal and the eclipses 
only partial, TSOO incorporated their spectroscopic light ratio as an external constraint 
on their light curve solution. In this way they derived a ratio of the radii of 1.05 ± 0.02 
(secondary/primary). The unconstrained solution by MOO, on the other hand, leads to a 
value of k around 1.5, similar to that obtained originally by DW35 also without the use of 
an external constraint. Such a large value seems somewhat unusual for main-sequence stars 
with very similar masses (M sec / 'M pi . im = 1.0279 ± 0.0031; TSOO). We note also that the 
maximum amplitude of the minima from the MOO solution, setting i = 90°, is ~0.45 mag 
(see Fig. 1), whereas their empirical fit to the measured amplitudes, mentioned earlier, 
suggests peak values considerably higher than this. Several of the measurements of the 
depth of the eclipses presented by the same authors (their Table 10 and Table 11) also reach 
larger values (as large as 0.76 in one case) between 1900 and 1910. Similar indications are 
seen in the work by Lehmann (1991). Admittedly some of these measurements suggesting 
deeper minima have rather large errors due to the photographic nature of the material, 
and two of them by MOO at about the epoch of maximum amplitude actually have much 
lower values, as do two independent measurements by Lehmann. Further light may be 
shed on this issue by examining the individual brightness estimates for SS Lac from the 
Harvard plates as re-measured by Schiller et al. (1996) (see also MOO). These measurements 
were kindly provided to the author by E. F. Milone. Fig. 2 shows a sequence of those 
measurements from 12 Harvard plates taken on the night of 1902 May 30 (JD 2,415,900) 
during a primary eclipse. The brightness estimates have been referred, for convenience, to 
the mean light level of the system outside of eclipse. Comparison with the same curves 
shown in Fig. 1 gives the clear impression that the primary minimum was indeed deeper 
than 0.45 mag, and that the data are considerably better represented by the TSOO curve. 
The solution by TSOO for i = 90° indicates depths of about 0.7 mag, which is close to what 
is expected for component stars of roughly equal brightness. 

The adoption of a representative light curve solution has a rather significant impact 



2 The star designations we adopt here follow the photometric convention of referring to the star eclipsed 
at phase 0.0 (Min I) as the "primary" . For SS Lac this happens to be the less massive and smaller star. The 
spectroscopic convention (e.g., TM98, TSOO) has the names reversed. 
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on the results described below. Although the TSOO solution is certainly not exempt from 
criticism, based on the arguments presented above we have chosen for the present analysis 
to adopt that solution as a more consistent representation of the physical properties of the 
stars in SS Lac, pending improvements that may become possible as we discuss later in §6. 
Using EBOP as described earlier we have determined the functions Dj(i) and D//(i) for the 
depths of the primary and secondary eclipse, respectively, which do not differ greatly. They 
are represented graphically in Fig. 3. 

2.2. Time variation of the inclination angle, i(t) 

The simplest approximation, and the one usually adopted, is that of a uniform variation 
of the inclination angle with time. As a convenient expression we adopt 

i(t) = 90° -di/dt (t-tgo) , (1) 

where tgo is the epoch at which eclipses are central (i = 90°). In §5 we describe a more 
realistic representation but we show also that, given the quality of the data at our disposal, 
the linear expression is a sufficiently accurate description in the time interval covered by the 
present data. For other purposes such as large extrapolations or determining the boundaries 
of the eclipse season of SS Lac, a linear change in i is inadequate (see §6). 

Several estimates of the slope di/dt have been made in recent years based on a variety 
of assumptions: 

1. Lehmann (1991) used measurements of the eclipse amplitudes based on photographic 
plates from the Sonneberg Observatory in a procedure similar to what we described above, 
and obtained di/dt = 0.18 ± 0.02 deg yr _1 . He also estimated the eclipses to have been 
central in 1911 ± 3. Few of the details are given, however, and the original data are not 
available except in graphical form. 

2. TM98 assumed typical radii for the stars (R = 2.25 Rq) based on an adopted 
spectral classification A2 V. From the linear separation between the components derived 
from their spectroscopic orbital solution they estimated the inclination angle at the end 
of the eclipsing season to be 83°. The authors assumed that the eclipses were central in 
1900-1915, and also that they disappeared at about 1960 based on Lehmann's (1991) work. 
From this they arrived at di/dt = 0.13 ± 0.01 deg yr _1 . Extrapolating to the mean epoch of 
their radial velocity measurements they estimated i(1998) ~ 78°. 

3. MOO estimated the inclination angle at the mean epoch of the velocity measurements 
by TM98 to be i(1998) = 76?5, through a procedure involving iterations between separate 



- 9- 



spectroscopic and light-curve solutions using the photometric data by DW35. Adopting an 
epoch of 1912 for central eclipses, they then obtained di/dt = 0.157 deg yr _1 . 

4. TS00 combined the photographic measurements by DW35 and Wachmann (1936) 
to obtain a light curve solution giving % = 87?6 ± 0?2, which they assigned to a mean epoch 
of 1912 ± 10. Based on this analysis and on their spectroscopic orbit, they estimated the 
minimum inclination angle for eclipses to occur to be i m { n = 81?6, and adopted 1951 as the 
epoch of cessation of eclipses. From this they derived di/dt = 0.151q;q| deg yr _1 . 

Although all these estimates of di/dt turn out to be rather similar, they are not 
independent and the assumptions made are quite varied and sometimes inconsistent, mostly 
due to limitations in the observational material. In the sections that follow we describe our 
own efforts based on observations that have not been used before for this purpose. 

3. The data 

Measurements of the depth of the eclipses in SS Lac were compiled by MOO from 
several sources: their own re-measurement of many of the original Harvard plates used by 
DW35; photographic observations published by Wachmann (1936), Nekrasova (1938), and 
Mossakovskaya (1993); visual determinations by Kordylewski et al. (1961); and more recent 
photoelectric measurements obtained by them at the Rothney Astrophysical Observatory 
and also by Mossakovskaya (1993), after the eclipses stopped. A total of 43 estimates of 
D were collected for the primary minimum, and 27 for the secondary minimum, over a 
period of nearly 100 years. Some of them are averages determined from several individual 
brightness measurements. In addition to these, Lehmann (1991) reported estimates of the 
amplitude of the eclipses from Sonneberg Observatory plates, but they were published 
only in graphical form. Further photographic observations of SS Lac were made by 
Tashpulatov (1965) from 1937 to 1955. However, these data have been called into question 
by Mossakovskaya (1993), and a re-examination of some of the original plate material 
by other investigators has cast similar doubts (see MOO). Because of this, and because 
they contain relatively few measurements at low light levels, we have not considered the 
Tashpulatov observations here. With regard to the other measurements, for the purpose 
of this paper we ignore the minor difference between the visual and photographic bands, 
following MOO. 

The crucial ingredient in the selection of brightness measurements for 
estimating the amplitude of the eclipses is the ephemeris used to establish 
which observations are near an eclipse. MOO adopted the ephemeris by TM98, 
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Min I = 2,450,716.32(±0.15) + 14.41638(±0.00010) • E, which was the most recent at 
the time. A much more accurate determination is now available from TS00, based on 
all published times of eclipse as well as their new radial-velocity measurements, and 
accounting also for light travel time (due to the presence of the third star) and apsidal 
motion. The sidereal period (P = 14.4161471 ± 0.0000089 days) is significantly different 
from the previous value, and one order of magnitude more precise. Over the time interval 
in which the eclipse amplitude measurements of MOO are non-zero (~1890-1940, or 50 yr), 
the difference in the periods alone accumulates to a phase shift of 0.02, almost exactly half 
of the width of an eclipse. In extreme measurement regarded by MOO to be near 

mid-eclipse may actually fall beyond the limits of eclipse, leading to an underestimate of 
depth of the corresponding minimum. This is especially true given that MOO adopted the 
conservative approach (precisely to allow for possible errors in the TM98 ephemeris) of 
including all brightness measurements they determined to be within ±0.04 in phase from 
the minima, which is actually twice the true half-width of the eclipses. 

An illustration of this effect can be seen in Fig. 4. In the top panel we show 17 
visual brightness estimates obtained near primary eclipse by Kordylewski et al. (1961) 
on the night of 1930 November 28 (filled symbols), phased with the ephemeris by TM98. 
This is the longest published sequence of measurements made on a single night during an 
eclipse that shows a clear change in brightness. The light curve displayed for comparison 
purposes is from the solution adopted in §2.1, with a maximum depth appropriate for the 
date of these observations. Also shown with a different symbol are all measurements by 
Wachmann (1936) that are near primary eclipse. These observations are confined to a 
2-yr interval (mean epoch ~1931.8), as mentioned in §2.1, that happens to be very close 
to the Kordylewski date. The phase shift is fairly obvious for both data sets. When the 
same measurements are folded with the ephemeris by TS00, on the other hand, the fit is 
considerably better (lower panel). We point out, incidentally, that both the Kordylewski 
observations and the Wachmann observations are completely independent of the ephemeris 
by TS00 since they were not used by the authors in its derivation, and thus they serve as 
an independent check on its accuracy^. 



3 The careful reader may have noticed that the Harvard observations displayed earlier in Fig. 2 
(JD 2,415,900, ^1902.4) were folded with the TM98 ephemeris, yet they show good phase agreement with 
the predicted light curves, seemingly conflicting with the evidence above. As it turns out, if the phasing is 
done instead with the TS00 ephemeris the result is essentially the same because both prescriptions made 
use of the main epoch of primary minimum from DW35 in their determination, which happens to be the 
very night of these measurements. In the case of TM98 this epoch is imposed exactly, and for TS00 it enters 
more weakly as part of a least squares adjustment. The main difference will then be simply a phase shift, 
equal to the residual of that epoch of primary minimum from the TS00 fit (see their Table 6). 
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Additional phase errors may be incurred because of uncertainties in the epoch, or even 
in the geometric elements e (eccentricity) and w, which essentially determine the location 
of the secondary minimum in relation to the primary through the quantity ecoscu. Fig. 1 
shows some hint of this effect (lower panel). 

As pointed out by MOO, it is to be expected that some measurements will fall on the 
branches of the minima, and therefore they will underestimate the true amplitude. To assess 
the importance of this effect they plotted the measured eclipse depths as a function of phase 
(their Fig. 6). From the fact that not all low amplitude values showed large displacements 
from the centers of the eclipses they concluded that the effect was unimportant. If the 
phases are computed instead with the new ephemeris, the general impression from such 
a plot is rather different (see Fig. 5): nearly all low amplitude values are quite far from 
mid-eclipse, and several are even beyond the true limits of the eclipse (c.f. Fig. 1). This 
suggests that some of the amplitudes may indeed be underestimated, and again that the 
contamination from measurements outside of eclipse cannot be ignored. 

This prompted us to revise the compilation of eclipse depths by MOO by re-examining 
all the original data using the new ephemeris. As a result, our list of amplitudes for 
the primary and for the secondary eclipse is somewhat different: some estimates were 
removed because they lie outside of eclipse, and new ones were added. The 591 Harvard 
measurements used by MOO were made available to us by the senior author, and in addition 
we considered the measurements from Wachmann (1936), Nekrasova (1938), Kordylewski 
et al. (1961), and Mossakovskaya (1993). 

For each source we first computed the average magnitude outside of eclipse. Each 
measurement during an eclipse was then compared to this average level to derive the 
amplitude. In the case of the Harvard data three observations during primary minimum 
and one during secondary minimum were much brighter than the average level outside of 
eclipse (by up to 0.75 mag), and were rejected. On nights with several brightness estimates 
during eclipse we averaged the data as well as the dates, except in cases when the sequence 
of measurements crossed over from the descending branch to the ascending branch. In 
such instances the average brightness will always underestimate the true amplitude at the 
average phase, and therefore we split the data at the center of the eclipse. The uncertainty 
we assigned to an amplitude estimate is the standard deviation of a single observation 
outside of eclipse, for each source. In cases where several measurements were averaged, 
the error was decreased by the square root of the number of observations. Finally, we 
included an estimate in 1991 derived from the Hipparcos measurements (ESA 1997), which 
was shown by TS00 to be very close to mid-secondary eclipse. Table 1 and Table 2 list 
all the measurements we considered for the primary (40) and the secondary eclipse (45), 
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respectively. 

Unfortunately the original dates of the measurements by Lehmann (1991) have not 
been published, and this prevents us from recomputing the phases to establish whether they 
can be incorporated into the analysis. We have therefore chosen not to use them here. 

4. The algorithm and results 

The procedure we followed consists of fitting the D(t) model (composed of Dj(i) 
and Djj{i) for the primary and secondary eclipse, and i(t), as described above) to the 
eclipse depth measurements, in order to determine i(t). Measured amplitudes for both the 
primary and secondary eclipses were used simultaneously. Because many of these estimates 
are not exactly centered at the minima, and therefore will always underestimate the true 
amplitude, they must first be corrected for this systematic effect. This was done as follows. 
For each date of observation we predicted an inclination angle using a preliminary estimate 
of i(t). With this we computed a synthetic light curve using EBOP, which therefore has an 
amplitude appropriate for the date of the measurement. Next we calculated the phase of 
the observation using the ephemeris, and we applied a correction to the measured depth 
equal to the difference in magnitude between a point on the curve with the same phase 
as the observation, and a point at the center of the eclipse. Amplitude measurements 
displaced in phase by 0.017 or more from the center of a primary or secondary minimum 
(half width = 0.02) were not considered, to guard against possible errors in the ephemeris 
(period or epoch, or perhaps even residual variations that have not been accounted for 
caused, for example, by the presence of additional stars in the system) and also because 
such a small drop in brightness is essentially within the errors of measurement of the visual 
and photographic techniques. We then carried out the fit for di/dt and tgo, and improved 
our preliminary estimate of i(t). The procedure is iterative, and at each step the corrections 
are slightly different. 

The eclipse ephemeris reported by TS00 was derived by those authors simultaneously 
with the apsidal motion and with the spectroscopic elements of the eclipsing binary and 
of the tertiary star (including light-time effects, which can amount to 1/10 of the eclipse 
half-width). Because the apsidal motion depends on the assumed inclination angle, and this 
angle is time-dependent, TS00 adopted their own determination of i(t) to account for this 
variation. Therefore, the phases we initially compute here from this ephemeris also depend 
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in principle on that determination of i(t)n In order to be self-consistent, for the present 
analysis we have recomputed the ephemeris at each iteration using the same methods and 
the same data employed by TSOO (times of eclipse, and radial velocities), starting from 
their estimate. Successive iterations use updated expressions for i(t), and the change in 
the phases calculated at each step was found to be very small. The apsidal motion and 
spectroscopic elements also did not change significantly. 

The procedure converged in 7 iterations, using standard non-linear least-squares 
techniques at each step to solve for di/dt and tgo (e.g., Press et al. 1992). In order to 
ensure that we converged at the absolute minimum of the x 2 surface, rather than at a local 
minimum, we checked the results using a genetic algorithm (see Charbonneau 1995), which 
explores all of parameter space. Fig. 6 displays our best fit to the amplitude measurements 
of the primary and secondary eclipse, in which the resulting time variation of the inclination 
angle is represented by i(t) = 90° — 0?1441(81)[t — 1905.6(1.4)]. The maximum amplitude 
(corresponding to central eclipses) according to this fit was therefore reached somewhere 
between 1904 and 1907. Table 1 and Table 2 include the phase displacement of each 
measurement, the final amplitude corrections we applied, and the corrected eclipse depths 
as used in the least squares adjustments and figures. 

5. A more realistic description of i{t) 

As pointed out earlier, a linear approximation to the time dependence of the inclination 
angle may be valid over relatively short time intervals (i.e., the duration of the observations 
in our case), but in reality the change in i is more complex. Perturbations induced by the 
third star on the inner pair cause the angular momentum vectors of both orbits to precess 
around the total angular momentum vector of the system, which is fixed in space. The 
phenomenon is known as the "regression of the nodes" effect (Soderhjelm 1975; Mazeh & 
Shaham 1976). One of the observable consequences is a periodic change in the inclination 



4 The TSOO determination of i(t) relied to some extent on observations by Tashpulatov (1965), through 
the estimated epoch of cessation of eclipses. Their ephemeris derivation indirectly incorporated this time 
dependence of i, although in practice it has little or no effect on the results. As mentioned in §3, the 
Tashpulatov observations are now considered to be unreliable. However, because of the procedure adopted 
here, any residual effect from those data on the present analysis disappears completely after the first iteration. 
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angle of the inner orbit that may be expressed conveniently for our case as 



cos % = cos I cos e — sin I sin e cos 



2tt ^ 

' ^90 + 75 {t — tgo j 

* node 



(2) 



in which / is the angle between the invariable plane of the system (which is perpendicular 
to the total angular momentum vector) and the plane of the sky. The angle e is measured 
between the orbital plane of the eclipsing binary and the invariable plane, and the argument 
of the cosine term on the far right represents the secular change in f2, the longitude of the 
ascending node of the binary orbit, with a period P no de- The angle / is constant, but in 
general e may also vary on long timescales, although this would be problematic to measure 
and is likely to be a small effect. For the purpose of this paper we will consider e not to vary. 
As before, tgo represents the epoch at which % reaches 90°. The angle f2go is the position 
angle of the ascending node at the same instant. From the condition of central eclipses Ogo 
is given by cosQgo = cot /cote, and thus the four unknowns in the new expression for i{t) 
are /, e, P node , and t 90 . 

Given the quality of the data available for SS Lac it is unrealistic to attempt to solve 
for all four parameters. However, the nodal period can be estimated from the theory of 
three-body interactions. TS00 gave an order-of-magnitude estimate of 500 yr. An improved 
value may be obtained from the expression derived by Soderhjelm (1975) (his eq.[27]), 
which depends not only on the masses of the three stars (and therefore on the unknown 
inclination angles of the inner and outer orbits), but also on the orbital periods (known), 
the outer eccentricity (known), and the relative angle between the orbits as well as the 
orientation of each orbit relative to the invariable plane of the system. Several useful 
constraints on some of these quantities were discussed by TS00. For example, from the fact 
that eclipses no longer occur the inclination angle of the binary at the present time must be 
smaller than % = 81?6, based on the size and separation of the components. In addition, the 
inclination of the plane of the outer orbit with respect to the line of sight is expected to be 
no smaller than 30° based on the spectroscopic solution, since smaller angles would imply 
that the third star would be massive enough and therefore bright enough that it would have 
been detected spectroscopically, if it is a main-sequence star. This, in turn, along with the 
measured apsidal motion, restricts the range of possible values for the relative inclination 
angle between the orbits (ei nner + e outC r) to be between about 24° and 30° (TS00). Though 
all these limits are helpful, the various angles involved are still not known exactly, so we 
have performed Monte Carlo simulations subject to all available constraints in order to 
determine the probability distribution of P no de- The result is shown in Fig. 7. Periods 
near 680 yr seem more likely, but the range of possible values extends from 600 to 700 yr. 
Soderhjelm (1982) has tested the accuracy of the analytical expression used above for P no d e 
by direct integration of the equations of motion in the three-body problem. For stars with 
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orbital and dynamical characteristics similar to those in SS Lac he concluded that P no de 
may be overestimated by that equation by up to 10% or so. We may therefore adopt as a 
representative value for SS Lac a nodal period of P no de = 600 yr. 

To arrive at the best fit to the observations we again used a genetic algorithm to 
search the parameter space, and found two solutions that are conjugates of each other 
(the result of a reflection of some of the angles around 180°). They represent physically 
equivalent scenarios that correspond to direct and retrograde motion of the binary on the 
plane of the sky. To determine which of these configurations holds for this system requires 
astrometric measurements, which are not available for SS Lac. In both cases the plane of 
the orbit of the binary is moving away from 90° at the present time, and the consequence 
for eclipses is the same. Mathematically the ambiguity is expressed by the inclination 
angle numerically increasing or decreasing from 90° during the 20 th century. Choosing 
the scenario in which i is decreasing (as was done implicitly in §4 and following the usual 
convention), we performed least-squares fits as in §4 to solve for J, e, and t 90 , using the 
data for both eclipses simultaneously. The corrections to the measured amplitudes were 
adopted from the final iteration in §4 (see Table 1 and Table 2). We obtained I = 84° ± 10°, 
e = 15° ± 3°, and tgo = 1905.8 ± 1.5, which are consistent with all available constraints. The 
angles / and e are significantly correlated, and this is reflected in their formal uncertainties. 
Fig. 8 shows this fit (solid lines) along with the corrected observations. For comparison we 
show also the solution from §4 using a linear model for i(t) (dashed lines), which is seen to 
be very similar to the new result and provides an equally good fit. The epoch of maximum 
eclipse amplitude is essentially the same in both solutions. 

6. Discussion 

Even though the fits from both models for i(t) are virtually indistinguishable given the 
errors in the observations, the corresponding predictions for the long-term behavior of the 
inclination angle and the implications for the eclipsing seasons of SS Lac are very different. 
According to the linear model the inclination angle completes a full revolution in 2500 yr, 
leading to eclipse seasons that repeat every 1250 yr. However, a secular change in i is not 
what is generally expected from three-body interactions. Instead, the inclination typically 
oscillates in a roughly sinusoidal fashion, the exact shape depending on the angles involved 
in eq.(2). The period of this modulation for SS Lac is about P no dc = 600 yr, as discussed 
above, which is the length of a full nodal cycle. The eclipse "seasons" (defined here as 
intervals when eclipses happen) can in general occur once or twice during the cycle. In the 
present case the indications are that there are two seasons per nodal cycle, each lasting 
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roughly 100 yr. This is shown in Fig. 9, where we display the depth of the eclipses (primary 
and secondary) over an interval of 1500 yr according to the model fit in §5, as well as the 
behavior of the inclination angle to the line of sight. Eclipses can take place only for values 
of % in the shaded area indicated in the lower panel (81?6 < % < 98?4; see TS00). The eclipse 
seasons recur in pairs, which combined last for about 350 yr. The observations during the 
20 th century (top two panels) are seen to cover the second season of one of such pairs, 
and no further eclipses are expected until shortly after the year 2200, if the present model 
is correct. For comparison we show also in the bottom panel the trend expected for the 
inclination angle according to the linear model (dashed line; §4), which is seen to match the 
slope of the more realistic model very well at the present time. 

The conjugate solution mentioned in §5 gives an identical pattern of eclipse seasons, 
but the function representing i(t) is a reflection about an axis at i = 90° of the curve drawn 
in the lower panel of Fig. 9. 

The end of an eclipse season depends in practice on the precision of the measurements, 
as argued by MOO. They concluded on the basis of their analysis of the amplitudes that, 
in the absence of photoelectric data, the eclipses effectively ended between 1937 and 1938. 
The slower decline resulting from the present fits favors a somewhat later date around 1945 
for the primary eclipses and ~1950 for the secondary eclipses, if we assume that amplitudes 
under 0.1 mag are undetectable by the visual and photographic techniques applied to 
SS Lac. Also, central eclipses occurred some 5 or 6 years earlier than in the MOO estimate 
according to the fits in §4 or §5. 

Though probably closer to the truth than the linear form of i(t), the more sophisticated 
"regression of the nodes" model is not without its weaknesses. The fit to this model 
is not as robust as the simpler approximation, partly because of the lack of amplitude 
measurements to the left of the maximum. Instead of two unknowns, three quantities must 
be solved for (after fixing P no de) out of the limited amount of information contained in the 
data. In addition, the errors given for /, e, and t 90 do not include uncertainties in the mean 
light curve adopted from TS00, in the value of -P no dc) or in the constraints used on e, which 
are based on the measured apsidal motion as reported by TS00. All these uncertainties 
are rather difficult to quantify. We cannot rule out that, given additional measurements 
of the depth of the minima, the detailed shape of the function depicted in Fig. 9 as well 
as the pattern of the eclipse seasons described above could change somewhat, particularly 
regarding the spacing of the two seasons in each nodal cycle. 

As it turns out, the inclination angle of the orbital plane of the binary extrapolated to 
the mean epoch of the spectroscopic observations by TS00 (~1998) remains ill-determined. 
From the linear model for i(t) we obtain 76?7, while the other model predicts 76?6, which 
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would seem to indicate good agreement. Although the formal errors in these quantities are 
roughly 1°, the true uncertainty is likely to be larger so that the absolute dimensions of 
the components can still not be computed reliably, especially considering also the possible 
systematic errors in the mean light curve (see below). 

The main limitation in the present analysis lies in the quantity and quality of the data, 
both for the amplitude measurements and the light curves. However, improvements can still 
be made in both areas. It is possible that useful information on the depths of the eclipses 
can still be extracted from the Harvard plates that were not re-measured by Schiller et al. 
(1996), as well as from the valuable plate material obtained at the Sonneberg Observatory 
and used by Lehmann (1991), but never published. The light curve solutions might also 
benefit from these new measurements, and this is important not only to improve our 
knowledge of the physical properties of the components in SS Lac, but for the determination 
of the depth functions Dj and D n as well (§2.1). 

The most difficult problem faced by all previous investigators attempting to fit the 
light curves of SS Lac has been the fact that the inclination angle of the binary changes 
with time. Inevitably this will affect the solution in subtle ways, and it is likely that 
several of the photometric elements are biased to some degree as a result of the averaging 
of the measurements during the minima, as the amplitude changes. A complication of a 
similar nature has been noted before by other investigators concerning the apsidal motion 
in eccentric binaries, particularly in systems with short apsidal periods. In this case, even if 
an accurate apparent period from a linear ephemeris is used for the primary, the effect is 
a distortion and widening of the secondary minimum because of the changing longitude of 
periastron. An example where this distortion is quite important is V477 Cyg. 

The solution to the problem of variable photometric elements would seem rather 
obvious, yet the common practice has been to circumvent the issue by dividing up the data 
into shorter time intervals, with the effect that data sets are poorer and the phase coverage 
more incomplete, as in SS Lac. If the change in the element (s) is linear to first order, as is 
often the case, it should not be difficult to correct for this directly in the computer program 
used to solve the light curve, allowing for the variation as a function of time. In the 
apsidal motion case, this has in fact been done for the system mentioned above, V477 Cyg, 
by Gimenez & Quintana (1992) using the computer program EBOP. They modified the 
appropriate subroutines so that the phase calculated for each observation accounts for the 
change in the longitude of periastron using a pre-determined value of du/dt. Similarly, 
either EBOP, the Wilson-Devinney code, or other programs in use could in principle be 
modified to account for di/dt in calculating the theoretical flux to be compared with each 
photometric measurement, thus referring all observations to a standard (arbitrary) epoch 
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with a fixed inclination angle. More ambitious investigators might be tempted to add di/dt 
as a variable to be determined from the solution along with the other light elements, but 
this may only be possible in the best observed systems. 

Alternatively, the individual photographic and visual observations for SS Lac may first 
be corrected along the lines described in §4 for the amplitude measurements, and then 
submitted to the light curve analysis. This would require no change in the current computer 
codes. The present analysis of the eclipse amplitude measurements could then be repeated 
once an improved light curve solution becomes available, and the procedure iterated until 
convergence. 

7. Final remarks 

Among the very few eclipsing binaries in which the amplitudes of the minima have been 
seen to change, to the point where the eclipses disappear completely due to the cumulative 
effect of the perturbations from a third star in the system, only a handful are as well 
documented as SS Lac. Despite obvious limitations in the observational material spanning 
more than a century, this has presented us with a rare opportunity to learn more about 
the configuration of the system in the light of current theoretical ideas about three-body 
interactions among stars. In particular, the eclipse amplitude measurements provide a 
valuable record of the change in the inclination angle that, with the possible exception of 
the succinct report by Lehmann (1991), has not been investigated in detail by applying the 
direct constraints offered by these observations. 

The three main contributions of this paper to the study of this phenomenon are (i) 
the use of a more accurate ephemeris (from TSOO) to select from the historical record 
the brightness measurements that were obtained during an eclipse; (ii) the application 
of corrections to the amplitude estimates to account for the systematically lower values 
typically measured on the branches of the minima; and (iii) the consideration of more 
physically realistic models than used before for the change in the eclipse depth, based on 
the known properties of the system (mean light curve solution) and expressions for the 
expected time variation of the inclination angle. For the latter we have explored both a 
linear approximation that is found to be quite satisfactory over the period covered by the 
observations, and a more rigorous prescription derived from the theory of the "regression of 
the nodes" effect in triple systems. 

Equally good fits are achieved with both expressions for i(t), but realistic long-term 
predictions on the future occurrence of eclipses in SS Lac can only be made from the second 
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approach. Our fits suggest that the eclipse seasons (duration ~100 yr) come in pairs that 
repeat after about 600 yr, and that we may expect the current lack of eclipses to last until 
approximately the beginning of the 23 rd century. 

We are grateful to the referee, Dr. E. F. Milone, for a number of helpful comments as 
well as for suggesting the use of, and generously providing the individual measurements 
of the Harvard patrol camera plates made by Schiller et al. (1996). This research has 
made use of the SIMBAD database, operated at CDS, Strasbourg, France, and of NASA's 
Astrophysics Data System Abstract Service. 
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Fig. 1. — Comparison between the light curves solutions for SS Lac by MOO (solid line) and 
TSOO (dashed line), with the inclination angle set to 90° to simulate central eclipses (see 
text). The lower panels show enlargements near the primary and secondary minima. The 
MOO curves were generated from the parameters listed by those authors and using a standard 
version of the Wilson-Devinney code, while the TSOO curve was computed using EBOP. 
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Fig. 2. — Individual brightness measurements for the night of 1902 May 30 (JD 2,415,900) 
from the Harvard plates (as re-measured by Schiller et al. 1996; see MOO), referred to the 
average light level outside of eclipse. The error bars adopted correspond to the uncertainty 
of a single measurement out of eclipse. The phasing of the data was done with the ephemeris 
by TM98. The light curve solutions displayed for comparison purposes are as in Fig. 1 (MOO: 
solid curve; TS00: dashed curve). The measurements indicate a primary eclipse amplitude 
larger than predicted by the MOO curve, and more consistent with the TS00 solution. 
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Fig. 3. — Calculated depth of the minima as a function of the inclination angle, for the 
primary (I) and secondary (II) eclipse. The light elements adopted are those of TSOO (see 
text). 
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Fig. 4. — Visual observations by Kordylewski et al. (1961) (filled circles) during primary 
eclipse on the night of 1930 November 28 (JD 2,426,309), and photographic measurements 
by Wachmann (1936) from a 2-yr interval with a mean epoch ~1931.8 (stellated symbols). 
The data are folded (a) with the ephemeris by TM98, which was adopted by MOO, and (b) 
with the new ephemeris by TSOO. Magnitude estimates are referred to the mean light level 
outside of eclipse. The typical uncertainty is shown in (b). Displayed for comparison in 
both panels are synthetic light curves based on the elements by TSOO, with an amplitude 
appropriate for the epoch of these observations. The phase discrepancy when using the 
TM98 ephemeris is apparent. 
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Fig. 5. — Eclipse depth measurements for SS Lac as tabulated by MOO, shown as a function 
of phase from the center of the primary (a) and secondary (b) minimum. The ephemeris used 
is that of TSOO. Most low amplitude levels are seen to be significantly displaced from the 
center of the eclipses, suggesting that they have been underestimated because they fall on the 
branches of the minima. Some of the points are seen to be beyond the limits of the eclipse, 
reflecting inaccuracies in the TM98 ephemeris used by MOO to select these measurements. 
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Fig. 6. — Fit to the measured eclipse amplitudes from Table 1 (photometric primary; panel 
a) and Table 2 (secondary, panel b), corrected to the center of each eclipse as described 
in the text. The fitting model shown here is a combination of the functions D(i) and i(t), 
with the latter assumed to be a linear function of time (eq.(l); see text). The symbols 
representing different sources for the measurements are as follows: crosses (Harvard; MOO), 
squares (Mossakovskaya 1993), asterisk (Nekrasova 1938; a single secondary measurement), 
filled circles (Wachmann 1936), open circles (Kordylewski et al. 1961), and open triangle 
(Hipparcos). 
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Fig. 7. — Probability distribution of the period of the nodal regression in SS Lac, from 
Monte Carlo simulations using eq.(27) by Soderhjelm (1975). All available constraints on 
the system have been used, as described in the text. 
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i(t) based on the "regression of the nodes" effect in triple systems (eq.(2); see text), rather 
than a uniform variation, as in Fig. 6. The latter fit is shown for reference (dashed lines). 
Symbols for the individual measurements are as in Fig. 6. 
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Fig. 9. — Long-term behavior of the eclipse amplitudes for the primary (a) and secondary 
(b) minima, resulting from the fit to the measurements based on the "regression of the 
nodes" effect. The cycle repeats with a period P no de = 600 yr (see text). Panel (c) displays 
the expected behavior of the inclination angle according to this model. The shaded area 
represents the range in which eclipses are possible, according to TSOO, and the dashed line 
shows the behavior based on the linear model for i{t). Both fits are seen to be essentially 
indistinguishable over the interval covered by the observations. 



Table 1. Primary eclipse amplitude measurements for SS Lac. 
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OQO 

UoU 


V) 
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23. 


Wachmann 


i 


2 


426,496 


499 


1931 


421 





162 


± 





059 


-0 


007 





063 





225 


24. 


Kordylewski 


3 


2 


426,525 


409 


1931 


500 





218 


± 





034 


-0 


001 





001 





219 


25. 


Wachmann 


1 


2 


426,568 


498 


1931 


618 





152 


± 





059 


-0 


012 





151 





303 


26. 


Wachmann 


1 


2 


426,626 


544 


1931 


777 





082 


± 





059 


+0 


014 





182 





264 


27. 


Kordylewski 


8 


2 


426,655 


260 


1931 


856 





186 


± 





021 


+0 


006 





048 





234 


28. 


Wachmann 


1 


2 


426,655 


264 


1931 


856 





372 


± 





059 


+0 


006 





048 





420 


29. 


Wachmann 


1 


2 


426,972 


300 


1932 


724 





302 


± 





059 


-0 


002 





005 





307 


30. 


Wachmann 


1 


2 


427,030 


215 


1932 


882 





132 


± 





059 


+0 


015 





190 





322 


31. 


Harvard 


2 


2 


427,361 


594 


1933 


789 





247 


± 





108 


+0 


002 





006 





253 


32. 


Kordylewski 


4 


2 


427,693 


327 


1934 


698 





026 


± 





030 


+0 


013 





142 





168 


33. 


Kordylewski 


2 


2 


428,039 


339 


1935 


645 





101 


± 





042 


+0 


015 





164 





265 


34. 


Kordylewski 


1 


2 


428,053 


469 


1935 


684 





061 


± 





059 


-0 


005 





026 





087 


35. 


Mossakovskaya .... 


1 


2 


428,082 


386 


1935 


763 





127 


± 





104 


+0 


001 





002 





129 


36. 


Kordylewski 


3 


2 


428,370 


532 


1936 


552 





167 


± 





034 


-0 


011 





096 





263 


37. 


Mossakovskaya 


4 


2 


432,767 


405 


1948 


590 





007 


± 





022 


-0 


013 





051 





058 


38. 


Mossakovskaya 


1 


2 


440,509 


203 


1969 


786 





247 


± 





104 


+0 


011 





000 





247 


39. 


Mossakovskaya .... 


1 


2 


441,619 


229 


1972 


825 





037 


± 





104 


+0 


010 





000 





037 


40. 


Mossakovskaya .... 


1 


2 


442,282 


410 


1974 


640 





077 


± 





104 


+0 


012 





000 





077 



a The amplitudes based on Harvard plates are from a re-measurement of that material by Schiller et al. (1996) 
(see MOO). 

b Phase displacement from the center of the eclipse, computed with the ephemeris by TS00. 
c Amplitude corrections (mag) computed as described in the text. 

d Corrected amplitudes (mag) as used in the least squares fits and figures in this paper. 
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Table 2. Secondary eclipse amplitude measurements for SS Lac. 





Source a 


N obs 




HJD 




Year 


Amplitude 


A<f> b 


AD C 


D 


d 

II 


1. 


Harvard 


1 


2 


412,679 


733 


1893 


593 





603 


± 





152 


-0 


003 





047 





650 


2. 


Harvard 


1 


2 


414,005 


891 


1897 


224 





064 


± 





152 


-0 


012 





420 





484 


3. 


Harvard 


1 


2 


414,034 


851 


1897 


303 





391 


± 





152 


-0 


004 





101 





492 


4. 


Harvard 


4 


2 


414,063 


808 


1897 


382 





386 


± 





074 


+0 


005 





157 





543 


5. 


Harvard 


1 


2 


414,106 


804 


1897 


500 





110 


± 





152 


-0 


012 





420 





530 


6. 


Harvard 


1 


2 


414,193 


637 


1897 


738 





205 


± 





152 


+0 


011 





408 





613 


7. 


Harvard 


3 


2 


414,207 


732 


1897 


776 





493 


± 





088 


-0 


011 





396 





889 


8. 


Harvard 


1 


2 


414,236 


678 


1897 


855 





576 


± 





152 


-0 


003 





066 





642 


9. 


Harvard 


1 


2 


414,265 


454 


1897 


934 





296 


± 





152 


-0 


007 





238 





534 


10. 


Harvard 


1 


2 


414,885 


586 


1899 


632 





606 


± 





152 


+0 


009 





353 





959 


11. 


Harvard 


3 


2 


415,029 


529 


1900 


026 





372 


± 





088 


-0 


006 





218 





590 


12. 


Harvard 


3 


2 


415,995 


534 


1902 


671 





463 


± 





088 


+0 


002 





070 





533 


13. 


Mossakovskaya . . . 


1 


2 


416,024 


357 


1902 


750 





587 


± 





104 


+0 


001 





026 





613 


14. 


Harvard 


2 


2 


416,370 


562 


1903 


698 





166 


± 





108 


+0 


016 





637 





803 


15. 


Harvard 


1 


2 


417,451 


546 


1906 


657 





633 


± 





152 


+0 


000 





000 





633 


16. 


Harvard 


1 


2 


423,621 


654 


1923 


550 





336 


± 





152 


-0 


003 





031 





367 


17. 


Nekrasova 


1 


2 


424,616 


350 


1926 


273 





383 


± 





193 


-0 


004 





046 





429 


18. 


Kordylewski 


1 


2 


425,409 


425 


1928 


445 





147 


± 





089 


+0 


008 





143 





290 


19. 


Harvard 


1 


2 


425,524 


647 


1928 


760 





193 


± 





152 


+0 


001 





004 





197 


20. 


Kordylewski 


2 


2 


426,332 


177 


1930 


971 





111 


± 





042 


+0 


016 





289 





400 


21. 


Kordylewski 


1 


2 


426,404 


242 


1931 


168 





021 


± 





059 


+0 


015 





274 





295 


22. 


Wachmann 


1 


2 


426,418 


661 


1931 


208 


-0 


038 


± 





060 


+0 


015 





274 





236 


23. 


Wachmann 


1 


2 


426,519 


494 


1931 


484 





152 


± 





059 


+0 


010 





172 





324 


24. 


Kordylewski 


3 


2 


426,634 


492 


1931 


799 





078 


± 





034 


-0 


013 





226 





304 


25. 


Harvard 


1 


2 


426,634 


615 


1931 


799 





055 


± 





152 


-0 


005 





051 





106 


26. 


Kordylewski 


1 


2 


426,649 


264 


1931 


839 





021 


± 





059 


+0 


011 





189 





210 


27. 


Wachmann 


2 


2 


426,649 


297 


1931 


839 





182 


± 





042 


+0 


014 





249 





431 


28. 


Wachmann 


1 


2 


426,692 


444 


1931 


957 





302 


± 





059 


+0 


007 





098 





400 


39. 


Kordylewski 


4 


2 


427,312 


364 


1933 


655 





148 


± 





030 


+0 


008 





113 





261 


30. 


Harvard 


1 


2 


427,600 


786 


1934 


444 


-0 


001 


± 





152 


+0 


015 





236 





235 


31. 


Kordylewski 


1 


2 


428,018 


480 


1935 


588 





021 


± 





059 


-0 


011 





158 





179 


32. 


TCordvlewski 


3 


2 


428 047 


417 


1935 


667 





114 


± 





034 


-0 


004 





028 


o 


142 


33. 


Kordylewski 


4 


2 


428,523 


262 


1936 


970 





311 


± 





030 


+0 


004 





028 





339 


34. 


Kordylewski 


2 


2 


428,552 


240 


1937 


049 





101 


± 





042 


+0 


014 





201 





302 


35. 


Kordylewski 


3 


2 


428,811 


362 


1937 


759 





081 


± 





034 


-0 


012 





161 





242 


36. 


Harvard 


4 


2 


428,811 


634 


1937 


759 





292 


± 





076 


+0 


007 





074 





366 


37. 


Mossakovskaya . . . 


1 


2 


429,849 


379 


1940 


601 


-0 


053 


± 





104 


-0 


008 





078 





025 


48. 


Mossakovskaya . . . 


1 


2 


433,006 


506 


1949 


244 





098 


± 





045 


-0 


010 





068 





166 


49. 


Mossakovskaya . . . 


1 


2 


433,900 


399 


1951 


692 


-0 


012 


± 





045 


-0 


004 





012 





000 


40. 


Mossakovskaya . . . 


1 


2 


440,503 


222 


1969 


769 





017 


± 





104 


+0 


010 





000 





017 


41. 


Mossakovskaya . . . 


1 


2 


441,238 


259 


1971 


782 





037 


± 





104 


-0 


003 





000 





037 


42. 


Mossakovskaya . . . 


1 


2 


442,305 


302 


1974 


703 


-0 


023 


± 





104 


+0 


014 





000 


-0 


023 


43. 


Mossakovskaya . . . 


1 


2 


442,636 


538 


1975 


610 





127 


± 





104 


-0 


009 





000 





127 


44. 


Mossakovskaya . . . 


2 


2 


447,826 


481 


1989 


819 





015 


± 





020 


-0 


001 





000 





015 


45. 


Hipparcos 


1 


2 


448,504 


165 


1991 


675 





001 


± 





025 


+0 


007 





000 





001 



a The amplitudes based on Harvard plates are from a re-measurement of that material by Schiller et al. (1996) 
(see MOO). 

b Phase displacement from the center of the eclipse, computed with the ephemeris by TS00. 
c Amplitude corrections (mag) computed as described in the text. 

d Corrected amplitudes (mag) as used in the least squares fits and figures in this paper. 



